Anomalous flux-flow dynamics in layered type-II superconductors at low temperatures 
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Low-temperature dissipation due to vortex motion in strongly anisotropic type-II superconductors 
with a moderate disorder (A 2 /Ef <C H/t <C A) is shown to be determined by the Zener-type 
transitions between the localized electronic states in the vortex core. Statistics of these levels is 
described by the random matrix ensemble of the class C defined recently by Atland and Zirnbauer 
flTfl , so the vortex motion leads naturally to the new example of a parametric statistics of energy 
levels. The flux-flow conductivity a xx is a bit lower than the quasiclassical one 0] and grows slowly 
with the increase of the electric field. 
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It is generally accepted that the flux-flow longitudinal 
(<j xx ) and Hall (a xy ) conductivities in the mixed state 
of type-II superconductors are given, respectively, by the 
Bardeen-Steven (BS) [GJ and Nozieres-Vincn Q relations 



Hc2 ~ i ^ 



(1) 



where a n is the normal-state conductivity, r is the elas- 
tic lifetime with respect to impurity scattering, and 
haj Ri A 2 /Ef <C A (it is supposed in (|]) that wo^ < !)• 
The meaning of this relation is very simple: it reflects the 
fact that the core region of a vortex (of the area ~ £ 2 ) 
may be considered (with respect to its electronic proper- 
ties) just as a normal metal. This idea is based on the 
existence of localized electronic levels within the core || , 
with energies constituting an equidistant set = phuj , 
where fi = i + 1/2 and i is integer. A microscopic theory 
of flux-flow conductivities, which explicitly takes into ac- 
count an existence of this level structure was developed 
long ago f§ (cf. also @, @], 0). This theory is based 
on the quasiclassical nonequilibrium diagram technique 
(see e.g. ||); its predictions for p xx and p xy basically 
coincide with the above simple picture as far as the su- 
perconductor is not in the so-called superclean limit (i.e. 
if the inverse electron- impurities scattering time r _1 is 
larger than the level separation loq). The underlying idea 
of this theory is that interaction between electrons lo- 
calized near the (moving) vortex and impurities can be 
treated as a kind of scattering problem in the continu- 
ous spectrum, i.e. similar to the impurity scattering in a 
normal state. 

In the present Letter we will show that in the case of 
strongly anisotropic (layered) superconductors the valid- 
ity of the above picture is limited to the situations when 
inelastic width T of the core electronic levels is compa- 
rable to or large than the level spacing loq. Note that 
r = Tjnj + T v , where T int is due to electron-electron and 
(mainly) electron-phonon interactions and grows with 
temperature, whereas r„ is due to nonstationarity of the 
impurity potential around moving vortex, and grows with 
its velocity v v ; therefore the condition r > wo is ful- 
filled at sufficiently large T or v v . In the opposite limit, 



r <C U)0i the core levels are essentially discrete; as a re- 
sult, all dissipation induced by the motion of a vortex 
is due to rare nonadiabatic (Zener) processes of electron 
excitations between the core levels and subsequent down- 
ward electron transitions with emission of phonons. As 
a result, the longitudinal flux- flow conductivity a xx be- 
comes reduced compared to its "quasiclassical" value (Q) , 
and dependent on the electric field: 



ecn e 

< ^^W T 1 - 





D 



where the characteristic electric field E* is very small, so 
in the domain of applicability of our theory E ^> E* 
(here p is the mass density of the crystal, s is some 
average sound velocity and c is the velocity of light). 
Note that according to Eq. (Q) conductivity grows with 
the increase of electric field (or current), which is the 
trend opposite to the one known near the classical vor- 
tex depinning transition ||. Our result (|^) is expected 
to be valid in the range of not very high current densi- 
ties j < jx = joikpl)" 1 ^ 2 , where I — vpr is the elastic 
mean- free-path, and jo ~ en e A/pp is the depairing cur- 
rent density. Note that the above result refers to the 
moderately clean case loq <C t _1 <C A/%, where the av- 
erage density of states within the core (v c (E)) is constant 
(at E > hu>o) as in normal metal. Therefore the effects 
we are discussing are of mesoscopic nature in the sense 
that they "feel" the structure of correlations of energy 
levels. Thus the range of existence of these effects disap- 
pears with the growth of UfI'- Jk cx (kpl)^ 1 ^ 2 ■ 

We are going to consider strongly anisotropic supercon- 
ductors with a high effective mass anisotropy m c /m a 3> 
1; it means very weak dispersion of the core level's 
energies as a function of k z , (E ll (k z ) — E a (0))/K ~ 

pujQ 7 — -f- «C Wo- It allows, as long as the low- lying levels 

with p <C m c /m a are relevant, to treat the core levels as 
purely discrete, which formally corresponds to a "pan- 
cake" vortex in 2D superconductor (one important point 
of our arguments where 3D nature of vortex is essential 
will be explained later). 
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At low temperatures any dissipation of energy related 
to the vortex motion is due to transitions between dis- 
crete electron states within the core. In the absence of 
electron-phonon interaction these levels have almost zero 
width (inelastic width due to electron-electron interac- 
tions for the system of discrete levels at low tempera- 
tures is even much weaker than the electron-phonon one, 
cf. Il0j ) , whereas their exact locations depend on the par- 
ticular realization of impurities in the region of the size 
~ £ around the vortex center r v . The most direct way to 
study the statistics of random energy levels is to employ 
powerful machinery of the Supersymmetric Sigma-Model 
jyj . Here, however, we prefer to use technically simpler 
(although less general) path, sufficient for the calcula- 
tion of dissipation rate in the case u>qt <C 1. This condi- 
tion ensures that the position of each level Ei is strongly 
modified with respect to the bare levels || of an ide- 
ally clean superconductor. The second important point 
is that the wavefunctions corresponding to all these levels 
are confined in the same area ~ £ 2 within the vortex core. 
Therefore it is quite natural to describe the distribution 
of levels in terms of an appropriate Random Matrix En- 
semble (RME) When vortex is moving, the realiza- 
tion of disorder in the core region is changing, so the core 
levels are moving up and down, presenting a new exam- 
ple of a "parametric level statistics" studied before with 
regard to disordered metallic grains and quantum dots 
(sec e.g. (jl3)). In our case the parameter X(t) govern- 
ing the evolution of levels is just the vortex coordinate, 
so one can try to make use of the relation [j^,[l6],[l3| be- 
tween the rate of energy dissipation per pancake vortex 
(within the applicability of the standard Kubo approach) 
and the mean-squared level "velocities" ((dEi/dX) 2 ): 

(dW/dt) Kubo = ^TrhC(0)(dX/dt) 2 , (3) 

where (3 = 1,2,4 for the system described by, re- 
spectively, orthogonal, unitary and symplectic Wigner- 
Dyson ensembles, and C(0) is the normalized dispersion 
of level "velocities", C(0) = ({dE t /dX) 2 )/{hcj Q ) 2 (here 
HllJq = (E i+ i—Ei)). Vortex solution breaks T-invariance, 
whereas leaves intact invariance with respect to spin ro- 
tations (here and below we neglect Zeeman spin splitting, 
which is weak as long as B <C H C 2). Therefore one could, 
at first sight, conclude that the relevant for our problem 
Wigner-Dyson ensemble is just the Unitary (U) one, with 
(3 = 2. In fact, the situation is a bit more complicated, 
because of the specific symmetry of the Bogolyubov-De 
Gennes equations determining the core levels: each of the 
positive-energy levels E^ = E > has its exact mirror 
counterpart EJ = —E. As a result, the RME relevant 
for our problem does not coincide exactly with Unitary 
or any other of the standard Wigner-Dyson ensembles. 
Fortunately, the general classification of RME related 
to the mixed superconductive-normal systems was devel- 
oped very recently by Atland and Zirnbauer (AZ) |t7| , 



who identified 4 different types of RME depending on 
the presence or absence of time- and spin-reversal sym- 
metries. The class C of the AZ classification just cor- 
responds to our problem at hand. The level statistics 
within the C class can be described via the multi-particle 
"wavef unction" of auxiliary "free fermions" in the same 
way as it was done for the standard Unitary ensemble 
]l4| ; the only difference is the condition of vanishing am- 
plitude of a single-particle "wavefunctions" in the origin 
of the energy axis E — 0. As a result, the level statistics 
within the C class coincides with that of the U class as 
far as highly excited states with Ei ^> Uujq are involved, 
but differs for low-lying states. 

It is very instructive to compare the dissipation rate 
obtained via (||) with the results of the quasiclassical ki- 
netic equation approach [jjj. As it will become clear soon, 
the applicability of the result (^) to the vortex dynam- 
ics is limited to the case when highly excited levels are 
important, so at this stage we can neglect difference be- 
tween C and U ensembles, and use (g) with (3 = 2. Com- 
paring the definition of the vortex damping coefficient rj 
{dW/dt = i]v v 2 ) with Eq. (§), one finds 

7rhC{0) = r\Kubo = irhn 2 D(u T), (4) 

where the second equality in (0) follows from the results 
of 0] in the limit lu^t <C 1, n2D = n e d is the areal density 
of electrons, d is the interlayer spacing, and we put sub- 
script "Kubo" to emphasize that the expression (0) has 
the same meaning and the domain of applicability as the 
Kubo-Greenwood formula. The damping coefficient r/ is 
proportional to the longitudinal flux-flow conductivity: 
&xx — r/ec/irhB. The parameter C(0) has the dimen- 
sion of inverse area; one can define characteristic length 
C = CT 1/2 (0) with the following meaning: when vortex 
moves over the length £, the characteristic displacements 
of the core levels inside it become of the order of the mean 
spacing ojq. Then Eq. (|4|) tells us that C = (r^D^o 7 ") -1 ^ 2 - 
Note that C decreases (i.e. sensitivity of the level posi- 
tions to the shift of a vortex increases) with the decrease 
of disorder, and becomes of the order of the Fermi wave- 
length at ujqT ~ 1, i.e. on the border of applicability of 
Eq. (Q). On the other hand, in the extremely disordered 
limit, UfI ~ 1, the length C would become of the order 
of the "clean" coherence length £o- 

There are two ways to understand the relation (j^): 
i) to consider open system with finite width of energy 
levels and use the standard Kubo-Greenwood approach 
]l3t , and ii) to work with strictly discrete levels but take 
into account the nonadiabatic transitions between time- 
dependent (due to variation of external parameter X) 
levels Ei(X), as it was done in The first ap- 

proach is clearly invalid in our case as far as the inelastic 
width r C uo- Following the second approach |l6| we 
find that Kubo-like expression (||) can be safely used if 
the characteristic frequency of level perturbations v v /C 
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is much higher than luq, i.e. at v v 3> vk = vo(kpl)^ 1 / 2 , 
where vo = A/ pp. On the other hand, at low vortex 
velocities v v <C vk and low temperatures the probabil- 
ity of Zener transitions between levels is exponentially 
low (~ cxp(— vk/vv)) when the intcrlcvcl spacing is of 
the order of its average value Tlloq. In this case dissi- 
pation occurs when in course of "level dynamics" the 
spacing between some pair of levels becomes very small, 
6E < Ez = Twjq{v v /vk) and Zener tunneling becomes 
probable. 

The crucial stage of the dissipation process is deter- 
mined by the non-adiabatic transitions of some electron 
from the highest negative-energy (filled) core level E-\ 
to the lowest positive-energy level Eq. The rate of these 
transitions can be calculated by the method pl| modi- 
fied for the C class of the AZ classification (details will 
be published elsewhere [H): 



oo oo 

Rq = v v JdA JdeN(e, A) exp 
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where N(e, A) = fgj [^^yj exp 
is the joint probability distribution for the energy e of 
the lowest positive level and its "velocity" A far from 
the "avoided crossing" region. Note that the rate Rq 
is just the factor 2 lower than the analogous result for 
the Unitary ensemble [ fTs] , which would also be valid in 
our problem if the transitions between highly exited lev- 
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els would be considered: R 



i|>l 



2i?n- In a similar 



problem treated in [|15| (where all transition rates were 
equal, Ri — R), the rate of energy absorption from the 
source was given by fruoR, which coincided exactly with 
the result for the Kubo regime, Eq. (||). This striking 
coincidence takes place, as was shown in [|l5 16 , only for 



the Unitary ensemble (ensembles with the level repulsion 
parameter f3 = 1,4 lead to velocity-dependent friction 
coefficient r\ cx v 13 / 2 ^ 1 ). 

In our case the situation is more complicated due to i- 
dependence of the transition rates Ri; as a result, the dis- 
sipation rate depends on the distribution function fi(t) 
for the electron's population of the i'th core level, which 
is determined by competition between Zener processes of 
excitation (the energy being absorbed from the source) 
and energy relaxation to the phonon "bath" caused by 
electron-phonon interaction. The rate of electron transi- 
tions between the core levels with energy difference Huj 
due to emission of 3D phonons can be estimated [jTJj as 
r(w) = Lo^fph, where 7 P h = v{uj T)Tiujln2D / 'ps 3 < 1 (here 
it is assumed that T < Tiuj, and v is the numerical co- 
efficient of order 1). Thus the kinetic equation for the 
distribution function ft(t) is: 



dfi(t) 
dt 



-Wb7pfc - fj] - EO'-O/jIl -fi}), (6) 



\j<i 
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which is a kind of discrete "diffusion equation" in the 
energy space. Dimensional estimates show that the 
characteristic width of the stationary solution ff is 

ichar ~ (i?oM)7p/i) 1/4 ~ ^Vy/vK^DEF/flUJ^T) 1 / 4 ; 

usually i c har is large since v v should be larger than 
the pinning-determined critical velocity v c , (estimate for 
Vc/vk will be given below), whereas the second factor in 
the above estimate is always very large. One can qualita- 
tively associate with i c har some "effective local temper- 
ature" T e ff(v v ) ~ hu>oi c har 3> Tiwo of the core-localized 
electrons. Due to large effective width of /?*, the en- 
ergy dissipation rate dW/dt is close to its value 2RqTiluq 
for the Unitary ensemble (all Ri — 2Rq), which is also 
the result in the Kubo regime, as mentioned above. In 
order to find correction r\ — r\Kubo we need to specify 
the rates Ri. It is very difficult to determine Ri for 
general i ~ 1 because of the: i) randomness of ener- 
gies where "avoided crossing" between levels i — 1 and i 
happens and ii) energy-dependence of the mean density 
of states for the C-ensemble. Thus we employ the sim- 
plest model with correct asymptotic behavior of rates: 
Ri^o — Roo — 2i?o- Within this model and at i c har 3> 1 
we get dW/dt = 2Ro(hcu + {E }{f - /_i)), which leads 
finally to the upper bound (cf. below) for the vortex 
damping coefficient r\ = (dW/dt) jv 2 : 



V 



< 1 



VKubo 



VK -1/4 



(7) 



Ri+l{fi - fi+l) + Ri{fi - fi-l) 



where v\ ~ 1. Eq. (^|) together with the relation between 
vortex velocity and electric field, v v = cE/B, leads to 
the announced in Eq. (j^) result for the conductivity o~ xx . 
Eq. (Q) is valid up to temperatures T < T e f f (v v ) ; at 
higher T the width of //' and damping coefficient are 
determined by temperature instead of vortex velocity. 
Although the final result for a xx is rather close to the 
one obtained within quasiclassical picture (where contin- 
uous spectrum of electron states is assumed) , the intrin- 
sic mechanism of dissipation is quite different. One can 
understand it in the following terms: the energy dissipa- 
tion rate is a product of the rate of inelastic transitions 
Rin by the characteristic amount of energy 5E trans- 
ferred at each transition; within our picture Ri„ is much 
lower than in the quasiclassic approach, whereas 5E is 
much larger by almost the same factor. This consider- 
ation suggests that Eq. (||) does not necessarily mean 
that the resistivity p xx — a xx as the quasiclassical re- 
sults ([!]) would tell us. The point is that Hall conduc- 
tivity <T x yi being nondissipative, does not vanish in the 
absence of nonadiabatic transitions and the Hall angle 
Oh = arctan((7 xl) / ! <j xx ) is not necessarily small even at 
LJ{)T <C 1. The calculation of a xy within our picture will 
be postponed for the future studies; here we just note 
that it can't be done within purely RMT approach, since 
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the latter neglects completely the correlations between 
energies of the core levels and their angular momenta. 
Such correlations (which survive under moderate disor- 
der) are irrelevant for the energy dissipation rate (which 
depends on the short-scale structure of level correlations 
only), but are crucial for the transverse nondissipative 
force acting on a vortex. 

We treat Eq. (Q) as an upper bound for the dissipation 
due to the following two reasons: i) we have used, when 
deriving Eq. (j?]), a model that overestimates transition 
rates R n =£o\ fi) we have employed purely classical mas- 
ter equation (^J) for the probability distribution function 
fi(t). It means that we neglected possible phase coher- 
ence of electron states on a time scale Rq 1 between two 
successive Zener transitions. It was shown in jl9| that in 
the opposite case of a strict phase coherence the system 
exhibits a kind of 1-dimensional localization (with the en- 
ergy E n playing the role of a "coordinate"); as a result, 
(dW/dt) st — 0. We expect these phase-coherence effects 
to be weak in the case of a 3D vortex composed from 
"pancakes", since the multiplicity of allowed electron 
states corresponding to each "level" (in fact, very narrow 
band) i means that the problem of energy-localization 
should be considered as quasi-lD, with a rather large 
number of "transverse channels" . Thus the correspond- 
ing "localization length" in the energy space Ie is ex- 
pected to be large compared to the characteristic width 
ichar of "classical" stationary distribution function, so 
the energy-localization effects are rather weak. However, 
these effects should be important for the case of a purely 
2D superconductive system. 

In the above calculations we have assumed that 
the vortex is moving with a constant velocity v v , 
which means that the driving current density j = 
v v erj /(nh cos 6 h) is higher than the pinning-induced crit- 
ical current density j c , so the corresponding "critical ve- 
locity" v c <C iv On the other hand, v v was assumed 
be small compared to vr — A/ftfc^/ 2 / 1 / 2 . Critical cur- 
rent density for the pinning of individual pancake vor- 
tex in layered superconductor on weak impurities can 
be estimated M as j c w jaiO.laimp/ld) 1 ^ 2 , where j is 
the depairing current density and <7i mp = (Irii) -1 is the 
electron-impurity scattering crossection (n^ is the impu- 
rity concentration). Using "pessimistic" estimate^ <C 1 
we get the compatibility condition for our theory in the 
form 5 • IQ^kpcjimp/ d <C {loot) 2 < 1. Thus the pro- 
posed scenario may be realized in the case of very weak 
impurities with a crossection small on atomic scale. Note 
however that the above condition may actually be less 
stringent due to the decrease of j c with the magnetic 
field already in the range B <C H c2 §, and also in the 
presence of a considerable Hall component of motion, 
tan Oh > 1- We believe that the necessary conditions 
may be fulfilled in the clean single crystals of layered su- 
perconductor NbSe2 where the main source of disorder 
is due to random substitution of small percentage of Nb 



by Ta, which is known to produce very weak impurity 
centers [^(J . Our results may also be relevant for layered 
HTSC, however in that case the role of d-wave pairing 
symmetry should be studied. 

In conclusions, we proposed new mechanism of flux- 
flow dissipation operative in layered superconductors 
at low temperatures. Dissipative conductivity a xx 
is found to be close to the classical results [^] 
but is slowly field-dependent. We suggest that, un- 
der the same conditions, the Hall conductivity a xy 
and the flux-flow resistivity p xx may differ consider- 
ably from their classical values. Useful discussions 
with B.L.Altshuler, S.Bhattacharya, V.B.Geshkenbein, 
N.B.Kopnin, V.E.Kravtsov, A.I.Larkin, V.M.Marikhin, 
V.P.Mineev and G.E.Volovik are gratefully acknowl- 
edged. This research was supported by the RFBR grant 
# 95-02-05720, DGA grant # 94-1189 (M.V.F.) and 
ISSEP grant # a96-568 (M.A.S.). 
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